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ABSTRACT 

We present hydrodynamical models for the Cassiopeia A (Cas A) supernova remnant and its ob- 
served jet / counter-jet system. We include the evolution of the progenitor's circumstellar medium, 
which is shaped by a slow red supergiant wind that is followed by a fast Wolf-Rayet (WR) wind. 

The main parameters of the simulations are the duration of the WR phase and the jet energy. We 
find that the jet is destroyed if the WR phase is sufficiently long and a massive circumstellar shell 
has formed. We therefore conclude that the WR phase must have been short (a few thousand yr), if 
present at all. Since the actual jet length of Cas A is not known we derive a lower limit for the jet 
energy, which is ~ 10 48 erg. We discuss the implications for the progenitor of Cas A and the nature 
of its explosion. 

Subject headings: hydrodynamics — ISM: jets and outflows — ISM: individual (Cassiopeia A) - 
supernova remnants 
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1. INTRODUCTION 

Over the last decade evidence has emerged that sug- 
gests that at least some core-collapse supernovae are in- 
trinsically non-spherically symmetric explosions. The ev- 
idence is strongest for supernovae of stars that have lost 
most of their outer ( hydrogen-rich) env elopes, i.e., the 
type Ib/c supernovae (|Wang et al.ll200ll ). For those ex- 
plosions, the inner layers are exposed early on, and asym- 
metries in the core more easily survive the interactions 
with the outer layers. This implies that departures from 
spherical symmetry originate from deep inside the explo- 
sion. 

A better understanding of the explosion geometries is 
needed to provide further insights into what powers core- 
collapse supernovae. In the canonical explosion model 
the explosion is driven by deposition of neutrino energy 
into the region just outside the proto-neutron star. How- 
ever, up to now, computer simulations of this core col- 
lapse do not self-consis tently predict supernova explo- 
sions (jJanka et al.l l2007). In those simulations the role of 
magnetic fiel ds and stellar rotation is us ually neglected. 
Accor ding to iKhokhlov et all (|1999D and IWheeler et all 
(2000), magnetic fields and rotation may play a crucial 
role in the explosion mechanism, and may lead to bipo- 
lar explosions. Additionally, other explosion mechanisms 
that are based on acoustic and hydrodynamic instabili- 
ties can result in a symmetric, albeit not necessarily bipo- 
lar ex plosions (e.g. lBlondin fc Shaw 2007: [Burrows et alJ 
l200l . However, if one considers the most energetic su- 
pernova explosions; those associated with the long du- 
ration gamma-ray bursts (LGRBs) , it is very likely that 
these explosions are trul y bipolar. The ass ociated super- 
novae are of type Ic, see iDella Valid §006) for a review. 

The engines that drive the explosions associated 
with LGRBs may or may not be related to those of 
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"normal" core-collapse supernovae. In the collapsar 
model, LGRBs are the resul t of black hole formation 
(jMacFadven fc Woo slcv 1999) and thus have a distinctly 
different engine from that producing normal supernovae. 
Alternative models that consider LGRBs to be pow- 
ered by_lii£hlj^nagnetic : _rapidly rotating, neutron stars 
(e.g. iThompson et al.|l2004D or by tra ns-relativistic blast 
waves m supernovae ~(|Tan et al.ll200lh . allow for a contin- 
uum of bipolarity and explosion energies. In those cases, 
the amount of rotation and magnetic field strength, and 
the line of sight, determine whether we observe a "nor- 
mal" supernova or one associated with a LGRB. This 
study aims to shed some light on the intermediate case 
of a supernova that shows distinct bipolarity, but is ener- 
getically in the range of regular supernovae and does not 
have relativistic ejecta. We try to provide some insight in 
the requirements on the energy in the asymmetric part of 
the supernova, i.e. it's "jets", and on the type of progen- 
itor that was responsible for the circumstellar medium 
(CSM) at time of explosion. 

Two likely examples of bipolar sup ernovae are known 
in the local neighborhood: SN1987A (jWang et alJl2002t ) 
and the supernova remnant (SNR) Cassiopeia A, the 
subject of this paper. The bipolarity of Cas A has 
only recently bee n establi shed from optica l (Feseni ^OOl : 



Fesen et alJl2006ft. X -rav (IVinkl l2004t iHwang et al.ll2004 : 
Laming et al.l 120061) and infrared (|Hines et al.ll2004l ) ob- 
servations. These observations show that, apart from the 
long known "jet" region in the northeast, a somewhat less 
prominent protrusion is located in the southwest (Fig. [I]). 
In Figure 1 we show in red the image of the jet as shown 
in IHwang et al.l (|2004[ K In order to show the jet in the 
context of the overall emission it is combined with im- 
ages in silicon (green) and radio (blue) . For details and a 



discus s ion on the jet and its abun dances we refer to Vinkl 
(|2004l ). IHwang et alJ p004). and lFesen et"afl (|2006f T~ 



The jets extend out to a radius of at least 3.8 pc, for 
a distance of 3.4 kpc (|Reed et al.lll995f) . The reason to 
believe that these jets are the results of a bipolar explo- 
sion, rathe r than being caused by a bipolar structure in 
the CSM (jBlondin et al.lll996h . is the distinct elemen- 
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Fig. 



Three color image showing the location of the 



jets (red, see IVinM [2004; Hwang ct al. 2004) with respect to 
the bright X-ray shell of ejecta (green, X-ray Si XIII emis- 
sion) and radio synchrotron emission (VLA archival data). The 
jet image is obtained by taking the ratio of Si XIII over 
Mg XI X-ray line emission. (Public dom ain image based on 
the 1 Ms Chandra observation of Cas A ([Hwang ct al] 120041) 



http://www.astro.uu.nl/~vinkj/casa_jet_si_radio.jpg ) 



tal abundance patterns in both jet regions, with the jet 
material coming from deeper layers inside the star. The 
X-ray and optical data indicate that the jet material is 
rich in oxygen burning products (Si, S, Ar, Ca), while it 
lacks carbon- and neon-burning products (O, Ne, Mg). 
This is the reason why the jet, of which the emissivity 
is weak compared to the rest of Cas A, stands out by 
taking the ratio of the Si/Mg line emission. 

Interestingly, the supernova that caused Cas A seems 
to have some shared characteristics with the supernovae 
associated with LGRBs and X-ray flashes. Cas A's pro- 
genitor probably had lost most of its hydrogen enve- 
lope, given the lack of hydrogen rich, optically identified, 
ejecta. Although we are not claiming that Cas A was an 
LGRB or even an X-ray flash, the two non-relativistic 
jets suggest that it may be related and there may be a 
continuum of bipolarity in the explosion of supernovae, 
thus providing a possible link between LGRBs and nor- 
mal sup ernovae. The to tal mass of shocked ejecta is 
2-4 Mq (|Vink et al.l Il996h . and the explosion energy is 
about a factor of t wo more than the canonic al explosion 
energy of 10 51 erg dHwang fc Lamind 12001 ). The total 
oxygen ejecta mass of 1-2 M Pl su ggests a main sequence 
mass of 18 - 22 M Q (jVinkj l2004f ) . These properties are 
reminiscent of the parameters derived for SN2006aj, the 
supernova associated with the X-ray flash XRF 060218 
(jMazzali et al.ll2006[ ). and simi lar to SN2003jd, th e one 
suggested to relate to a LGRB dValenti et aLlfeoOTT) . 

The large amount of swept up mass in Cas A and 
the dynamic properties of the blast wave suggest that 
the blast wave is currently moving through the high- 
density red supe rgiant (RSG) wind (jChevalier fc OishU 
l2003HVinkir2004D . However, the lack of H-rich ejecta sug- 
gests that Cas A exploded as a Wolf-Rayet (WR) star. 
Moreover, the presence of slow moving N-rich knots has 
been explained as originating from the hydrodynamical 
instabilities between the fast WR wind and the d ense, 
slow moving, RSG wind jG arcia- Segura et al1ll996[ ). 

In this paper we present hydrodynamical simulations 
of the jets in the context of the progenitor's mass loss his- 
tory, which we take to be a RSG phase, possibly followed 
by a WR phase. There are two main reasons for pursu- 
ing this problem. First of all, the energetics of the jets 



can be better estimated using a realistic mass loss his- 
tory in the hydro-simulations. Secondly, the survival of 
the jets depends strongly on the mass-loss history of the 
progenitor. Therefore, the jets in Cas A can be used as 
a diagnostic on both the properties of the bipolar explo- 
sion, and on the progenitor-shaped circumstellar medium 
(CSM) at the time of explosion. 

2. METHODS 

The simulations of the supernova explosion and the 
rogenitor win ds are done using the code ZEUS-3D v3 . 4 



19961). an exten d ed ve rsion of the ZEUS-2D 



Stone fc Norman] (jl992). This code solves the 



Clarkd 
code by 

ideal non-relativistic fluid equations in three dimen- 
sions by finite differencing on an Eulerian mesh. Von- 
Neumann Richtmyer artificial viscosity is used to deal 
with shocks. Energy losses by radiative cooling are cal- 
culated according to the co oling curve as described by 
iMacDonald fc Bailevl (|1981[). Photo-ionizat ion is used as 
implemented bv lGarcia-Segura et alj (|1999l ). where mat- 
ter within the Stromgren radius is fully ionized, whereas 
all the rest is considered neutral. 

In our case the code is set up in spherical (r6<j)) coordi- 
nates with the assumption of symmetry around the polar 
0-axis. Since the supernova remnant does not reach be- 
yond the radius where the RSG wind meets the main 
sequence bubble, the initial grid is filled with a RSG 
wind. This part of the simulation is done in ID in order 
to save on computational time. Once the RSG wind oc- 
cupies the whole grid, the resulting CSM is transferred 
onto a 2D grid with 900 cells in the radial direction and 
225 cells in the angular direction. The equidistant grid 
cells are distributed over a radial extent of 6 pc and an 
angular extent of 7r/4 measured from the pole, meaning 
an effective resolution of 2.1 x 10 16 cm by 0.2°. The ap- 
plied resolution appears to be adequate for our purposes: 
doubling the resolution did not affect the hydrodynam- 
ical evolution of the structure and development of the 
instabilities. 

Into this grid, we release a WR wind that lasts for a 
period of 5.5 x 10 11 s (~ 17, 000 yr). During the simula- 
tion of this evolutionary stage of the supernova progen- 
itor, every 5.5 x 10 9 s a snapshot is taken of the CSM. 
These snapshots, representing the CSM for different du- 
rations of the WR phase of the progenitor, serve as the 
initial conditions in which we subsequently add the su- 
pernova ejecta and follow the evolution of the SNR. In 
other words, the age of the WR phase at the time of 
the explosion is a variable in the different simulations of 
the SNR, depending on the snapshot we take of the pro- 
genitor evolution to serve as the initial condition. The 
maximum duration of the Wolf-Rayet phase is restricted 
by the high amount of mass in the Cas A remnant and 
the high density the blast wave is currently running into, 
indicating that the swept-up shell, hereafter called WR 
shell, is within the border of the remnant (i.e. 2.5 pc). 

The parameters for the RSG and WR winds are taken 
from van Veelen & Langer (in prep. 2008). The prop- 
erties of the RSG wind are determined by the amount 
of mass lost by the progenit or star and the typic al RSG 
life time as in the model bv lHirschi et all (|2004l ). This 
gives a mass loss rate and terminal wind velocity of 
M = 1.54 x 10~ 5 MQyr" 1 and v HSG = 4.7 km s -1 . 
The temperature of the RSG wind is set to the typi- 
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cal value for this type of wind: T = 10 3 K. The prop- 
erties of the WR wind: M = 9.7 x 1(T 6 M yr" 1 , 
v = 1.7 x 10 3 km s " 1 , and T = 10 4 K, we re calculated 
in accordance with iNugis fc Lamersl (|2000D . Note that 
qualitatively our results do not depend on the exact pa- 
rameter chosen here, as will be discussed in section 13.31 
Observational co nstraints for the shocks are taken from 
several authors (iDeLanev fc Rudnickl l2003t IVink et al.l 
U998tlMorse et al.ll2004 : iKoraleskv et al.lll998h to be the 
following: forward « 5000 km s -1 , rf or ward ~ 2.5 pc, 
and the reverse shock to: ^reverse ~ 3000 km s _1 , 
1.4 pc. 



' reverse 



The explosion para meters are take n from 
iHwang fc Lamind (|2003h and IVinkl (|2004[ ) to be 
E ei = 2 x 10 51 , M ej =2-4 M , where we adopt an 
ejecta mass of M e j = 2.5 M Q . The initial ejecta density 
profile covers 0.1 pc in radius, and consists of a constant- 
density core, with an envelope for which the density 
decreases as p oc r~ 9 , which is the ty pical density profile 
of the ejecta in explosion models (c.f. lTruelove fc McKed 
1999). In order to match the observationally determined 
ejecta mass and energy, we iteratively determine the 
value for the density in the ejecta core and the velocity 
at which the core ends and the powerlaw envelope 
begins. For the parameters adopted in our simulations, 
this happens for a central density of the ejecta core of 
1.1 x 10 _19 g cm~ 3 , and the powerlaw envelope begins 
where the velocity of the ejecta is 9.8 x 10 8 cm s _1 . 
The velocity linearly increases from zero at the core, to 
15,000 km s _1 at the outer part of the envelope. 

Since the details of the jet-forming mechanism are 
poorly understood, we try to explore the relevant pa- 
rameter space in density and velocity by adjusting them 
in the jet region and thus locally enhancing the kinetic 
energy. Since it is not to be expected that the pre- 
supernova density will be much higher in the jet region, 
we focus on the velocity, although we do explore a factor 
of two in density decrease and enhancement. The shapes 
of the density and velocity profiles are expected to be the 
same as in the rest of the ejecta, since these are deter- 
mined by the p ropagation of the explosi on through the 
progenitor star (jMatzner fc McKed fl999T ). 

3. RESULTS 

During the progenitor's evolution, the CSM is shaped 
by a RSG wind that may be followed by a WR wind. 
The WR wind sweeps up a shell, whose radius and mass 
depend on the duration of the WR phase. In Figure[2]thc 
result of a 2D simulation of this evolution is plotted in 
ID and 2D. The swept-up shell is irregular due to hydro- 
dynamic instabilities. Since both winds are steady, the 
shell moves with constant velocity and is not Rayleigh- 
Taylor unstable. However, because of the large velocity 
difference between the WR wind and the RSG wind, and 
because cooling is implemented in the code and effective 
in the shell, the shell becomes sufficiently thin to be sus- 
cepti ble to the Vishniac "thin shell" instability (|Vishniad 
1983). This instability arises when the shell is thin 
enough and a compression ratio of over ~ 21 is reached, 
which is the case in our simulations. At late times in the 
WR shell evolution, the shell becomes more and more 
fragmented. However, this only occurs beyond the point 
we consider here, which is limited by the requirement 
that the SNR is currently running into the RSG wind 



and the shell therefore must be within a radius of 2.5 pc. 
For a more detailed descript ion of the instabilities and 
when t hey occur, we refer tolGarcia-Segura fc Mac Low! 
(|1995ft : iGarria-Segura et all ((19961). We find that dou- 
bling the resolution in both r and 9, reaching an effec- 
tive resolution of 1.0 x 10 16 cm by 0.1°, does not make a 
qualitative difference in the appearance of the instabili- 
ties and the overall dynamics. 

Figure EJa shows an average over all angles, which 
smooths out the density and pressure jumps. The free 
streaming RSG and WR wind follow the density pro- 
files (^p — M / (4irr 2 v W i n d)j that depend on their respec- 
tive mass loss rate and velocity. In the swept up shell 
the pressure is much higher than outside the shocked 
region. The shell consists, farthest from the star, of 
shocked, compressed RSG wind material. Further in- 
ward, the contact discontinuity, characterized by a jump 
in the density but with a constant pressure, separates 
the shocked RSG wind from the shocked WR wind. Fig- 
ure^ shows the density of the CSM in 2D. The Vishniac 
thin shell instability causes irregularities in the swept up 
shell. The three panels in Figure correspond to the 
last three evolutionary stages that were plotted in Fig- 
ure [2]a. Some accumulation of mass in the shell near 
the polar axis is due to the boundary conditions. This 
however does not qualitatively influence the outcome of 
the subsequent supernova simulation, as we tested with 
a simulation of the winds in the r</>-plane. 

The CSM that is created in this manner provides the 
background in which we implement the asymmetric su- 
pernova ejecta. The main parameters that we vary are 
the duration of the WR phase of our background and the 
energy in the jet. In the following part we focus on the 
results of different jet energies and use the background 
of a plain RSG wind. 



3.1. Jet Evolution in a RSG Wind 

With the explosion parameters as given in section O 
we find for the equatorial forward and reverse shock the 
following properties: ^forward = 5430 km s _1 , rf orward = 

2.43 pc, and Ureverse = 3160 km S^ 1 , Reverse = 1-63 pc, 

for an age of the remnant of 330 yr. The jet is colli- 
mated by the high pressure cocoon that results from its 
bow shock. The initial opening angle is taken to be very 
small (half opening angle of about 1°). Variations of this 
opening angle have been tried and we find that as long 
as it is not too large (< 4°), the half opening angle of 
the jet after 330 years is a lways about 10° (si milar to 
what is observed in Cas A, lLaming et al.ll2006T ). In nu- 
merous low-resolution simulations we have explored the 
effect of different jet densities and velocities on the re- 
sulting jet length. These results are shown in Figure [3l 
Differences in jet length for the same jet energy derive 
from the relative contribution of mass and velocity. For 
a given energy, a higher density gives rise to a longer jet. 

In order to see how the length of the jet scales with 
its energy, we have derived analytical solutions for both 
a cone-shaped jet with constant opening angle and for 
a cylindrically shaped jet (Appendix |A|) . These give the 
following relations between the length of a jet and its 
energy, for a cone: 
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Fig. 2. — a) CSM density (solid lines) and pressure (dotted) due to the RSG and WR winds are shown, for different stages in the 
evolution. The black solid line shows only the free-streaming RSG wind. For longer WR phases the WR shell broadens and the density 
contrast increases. The shocks appear smooth due to averaging over the 7r/4-angle we used in the 2D simulation, b) CSM density evolution 
is plotted in 2D, which exemplifies the irregularities in the shell due to the Vishniac thin shell instability. The three panels show the CSM 
for different durations of the WR phase and correspond to the last three profiles in a) . The forward shock is marked by the sudden increase 
in density and corresponds to the outer boundary of the shell, at a radius of approximately 0.35, 0.5, and 0.9 pc. The contact discontinuity 
corresponds to the sudden drop in density, and the inner boundary of the thin shell. The wind-termination shock is located at around 0.15, 
0.25 and 0.45 in the three panels respectively. 



i?= 1.15 



Et 2 
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where 6j is the opening angle, E the jet energy, R the 
length of the jet at time t, and A is a density normaliza- 
tion constant. For a cylinder the relation is: 



data from simulations 

location of outermost observed knot 



R = 



I Et 2 R 
C^Airr 2 ' 



(2) 



with r the radius of the cone and C 7 = (47)/((7 — 1)(7 + 
l) 2 ) = 45/32 for a polytropic index of 7 = 5/3. Neither 
of these relations gives a satisfactory fit to the simulated 
jet length. Clearly, the approximations used are too sim- 
ple to encompass all the relevant physics. For example, 
the assumptions arc based on a point explosion and ne- 
glect the presence of a reverse shock, lateral spreading 
of the jet and interaction with the rest of the supernova, 
which are also important for the propagation of the jet. 
However, it gives us an estimate of how the propagation 
of the jet may scale with energy and the order of magni- 
tude we can expect. 

In the "jet" -region, optical knots have been observed 
out to a dis tance of 3.8 pc fro m the center of the SNR 
(|Fesenll2001[ ). As mentioned bv lLaming et alj (|2006t ) this 
may not be the actual tip of the jet, which may be invis- 
ible due to cooling. The actual blast wave may be out- 
side of the field of view of available X-ray observations. 
The same authors estimate the tip to be at a distance of 
5.66 pc, with an equivalent isotropic energy for the jet of 
2.3 x 10 52 erg. Our findings concur with these estimates. 
For a jet length of 3.8 — 6 pc, our simulations show that 
the energy of the jet should be in between 6.0 x 10 47 erg 
and 5.0 x 10 48 erg, which is equivalent to 8.0 x 10 51 erg 
and 6.6 x 10 52 erg isotropic energy. 

3.2. The Evolution of a Jet in the Presence of a WR 

Shell 

The presence of a WR shell severely alters the mor- 
phology of the SNR and limits the survival of a jet. From 
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Fig. 3. — Data from simulations (black dots), representing the 
length of the jet for different jet energies. We have varied the initial 
opening angle, density and, mostly, the velocity, thus sometimes 
getting different jet lengths for similar energies. If the mass in the 
jet is relatively high, the jet protrudes further. The horizontal line 
indicates where the outermost observed ejecta are located. 

low-resolution simulations we find that for WR shells fur- 
ther out than ~ 1 pc the forward shock of the SNR ex- 
tends beyond the observed 2.5 pc. We therefore focus 
our higher resolution simulations on the earliest stages 
of the WR phase. 

In Figure HJa the logarithm of the density is plotted in 
2D. The three panels show the evolution of the same, ax- 
isymmetric, supernova ejecta into different CSMs. The 
CSMs vary only in the duration of the WR phase. The 
left plot shows the density after explosion into a plain 
RSG wind. In the middle and right plots, the progenitor 
star model did include a WR phase with a duration of 
respectively 2265 and 3480 years. The average density 
profiles of the different CSMs were shown in Figure EJa. 
We find that the jet does not protrude through the shell 
if the CSM contains a thick, high density, and therefore 
massive shell. This is the case in the right panel of Fig- 
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ure [4)3. The jet can punch through a less massive shell, 
but it will be shorter than in the case where the CSM 
was shaped solely by a RSG wind. The forward shock of 
the main remnant initially is perturbed when it collides 
with the circumstellar shell, but the shock smoothens 
within about 100 years due to the high sound speed (in 
the earliest phase of the remnant c s w 10 8 — 10 9 cm s ). 

The thermal emission due to bremsstrahlung is plot- 
ted in Figure [3^. Note that the thermal bremsstrahlung 
emissivity is only an approximation. In reality the X-ray 
emissivity is determined not only by the temperature and 
the density, but also by the composition of ejecta and 
CSM, equilibration of electron- and ion-temperatures, 
and non-equilibrium ionization, all of which are not taken 
into account. However, Figure|4^ shows that the jet emis- 
sivity is small compared to that of the high density parts 
of the shell, as is indeed what we see in Cas A. In the left 
and middle panel, most of the thermal emission of the jet 
is concentrated in the lower first half of the jet. Likely 
the outer part of the jet would not be bright enough 
to see in the observations. Therefore, the observations 
will probably not show the maximal extent of the jet, 
but rather only a part of it. The shock at the tip of 
the jet may be seen in non-thermal X-ray emission, but 
currently there are no detailed X-ray observations of this 
region. Because of the high temperature near the forward 
shock, the high density material in the outer part of the 
remnant stands out compared to the pure-density plot. 
What appears as a bright shell in the observations and 
is often cited as the reverse shock may therefore rather 
be closer to the contact discontinuity, around the higher- 
density fingers. Tracking of the different fluids would be 
required to determine exactly where the ejecta and the 
CSM meet. For a short WR phase (middle and right pan- 
els) the dumpiness appears that could be responsible for 
the observed bright nitrogen-rich knots in Cas A. This 
is in contrast with what happens for a smooth circum- 
stellar medium shaped solely by a RSG wind, as shown 
in the left panels. In fact, the contact discontinuity in 
this case is Rayleigh- Taylor unstable, which, because of 
limited resolution, does not show up unless we initially 
add 1% density perturbations in the RSG wind. 

3.3. Shell Density and Jet Survival 

In the previous section, we showed that the jet remains 
present in the SNR only if the progenitor has a very 
short WR phase, while it is stalled for a progenitor with 
a longer WR phase. To understand the physics behind 
this, we compare the energy of the jet with the energy 
that is required to accelerate the part of the shell within 
the opening angle of the jet to typical post-shock veloc- 
ities. The jet does not significantly slow down before it 
reaches the shell, so for the post-shock velocity we take 
3/4 of the initial jet velocity Vj. The mass in the shell is 
calculated by evaluating the density and the volume of 
each grid cell that lies within the forward and the reverse 
shock of the shell and within the opening angle of the jet 
0j, and correct it for the solid angle. The energy that 
is needed to accelerate that portion of the shell to 3uj/4 
should then be: 

E acc = ~M^(Sl i )(~y?\ . (3) 



Here M s h e ii(^j) represents the mass of the shell within 
the solid opening angle of the jet. 

The "required acceleration energy" is plotted in Fig- 
ure[5]for three simulations of the earliest phase of the WR 
shell, and compared to the energy in the jet. It appears 
that this gives a good measure for determining whether 
the jet will remain collimated after the encounter with 
the shell, or not. When the energy in the jet is higher 
than the energy that is required to accelerate the portion 
of the shell, it protrudes, whereas when the acceleration 
energy needed is much higher, the jet will disperse into 
the rest of the remnant. It is not a perfect indication; 
the jet can survive when the energy needed to accelerate 
the shell is up to a few times higher, although the length 
of the protruding jet will be limited. 

The dotted portions of the curves indicate the situation 
where the jet does not protrude, whereas the solid and 
dashed curves indicate the region where the jet is larger 
than ~ 3.8 pc. In order to rule out a coincidental correla- 
tion, we did three different simulations of the shell in the 
WR phase. We varied the initial density perturbations 
in the different runs. As a result, the density fluctuations 
in the shell, and hence the mass accumulation, differ for 
the three cases. On average, the mass in the shell should 
be approximately equal to the mass of the RSG wind 
out to a radius equal to the outer edge of the shell; since 
the velocity of the WR wind is much higher than that of 
the RSG wind, the density of the WR wind is lower and 
it's mass contribution negligable compared to the mass 
of the RSG that is swept up by it. Therefore, on average, 
the total energy required to accelerate the shell increases 
with the time as the WR shell is allowed to develop. 
The dotted parts of the curves indicate the time where 
the jet does not have enough energy to break through 
the shell. The curve that starts as a solid line represents 
the WR shell simulation that was plotted in Figure [2] and 
was used for the SNR/jet simulations that were depicted 
in Figure [5] The horizontal dashed line indicates the 
energy of the jet used in the simulations. The vertical 
patch of line indicates the analytical solution where the 
mass, were it to accumulate homogeneously, i.e. without 
instabilities and equal for each unit angle, multiplied by 
i(|«j) 2 equals the jet-energy. More explicitely, where: 



£ jet = -M r8g (iWi,fij) ( (4) 
= ^~a • ( 5 ) 

47TWRSG 

We can see that for a WR phase of about 2000 years, a 
jet energy of less than 10 49 erg (~ 1.3 x 10 53 erg isotropic 
energy), is not sufficient. The exact time varies with each 
simulation, due to small random differences in the devel- 
opment of the WR shell. However, in all three cases, the 
turning point lies close to the analytical value of about 
1800 yr. The result of the jet-development we have seen 
in Figure [4] Indeed the jet protrudes in the case of a 
very short WR phase progenitor, and is stalled when the 
progenitor's WR phase is longer than about 2000 yr. 

Note that the critical duration of the WR phase for 
which a shell forms that is massive enough to block the 
jet, depends on the exact properties of the progenitor 
winds. Although in our model the properties of the RSG 
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a) 1VK phB.se= 0, 



8B5, 



j I WR phaae= 0. 




Fig. 4. — a) Density and b) approximate X-ray bremsstrahlung emissivity of the supernova remnant at a time of 330 yr after explosion. 
The left panels of the density and X-ray figures show the remnant that results from the evolution of asymmetric ejecta in a CSM without 
a Wolf-Rayet phase, i.e. a pure RSG wind. The middle plots show the case where the Wolf-Rayet phase has lasted 2265 yr, and the right 
panels show the remnant of an explosion into a CSM with a Wolf-Rayet phase of 3480 yr. In all of these cases the density of the jet is 
equal to the density in the rest of the SN, and the velocity is enhanced with a factor 6, giving a jet energy of 5.5 X 10 48 erg with a maximal 
velocity of 90 Mm s — 1 . The forward shock of the remnant in all three cases is located at a distance of ~ 2.4 pc. The reverse shock, however, 
is located increasingly farther inward for longer WR phases, due to the extra pressure created by collision with the shell. In the left panel 
the reverse shock is at a distance of ~ 1.6 pc, whereas for the middle and right panels it is located at around 1.0 and 0.8 pc. The contact 
discontinuity of the main remnant for the left-most plots is located at a distance of 1.8 pc, right in front of the high density shell. The 
contact discontinuity in the middle and right-hand panels is rugged due to the presence of the shell, but is on average located near the 
high-density fingers. For calculating the X-ray emissivity, the electron temperature is set to be equal to the plasma temperature. Note that 
the emissivity is not integrated over the line of sight but represents a slice through the meridional plane. 

wind are well constraint by the observations, the proper- r 

ties of the WR wind are more uncertain. The physically 
important factor is the mass contained in the swept up 
shell. For a higher mass loss rate or higher velocity of the 
WR wind, the critical duration of the WR phase may be 
found to be somewhat shorter. A lower mass loss rate or 
lower velocity will result in later accumulation of the re- 
quired mass and induce a longer critical duration of the 
WR phase. We simulated jets with energies in excess 
of 10 49 erg, but find that when the jet energy becomes 
comparable to the explosion energy, the jet and the SNR 
shell broaden and no longer resemble the morphology of 
Cas A. 

In the simulations, the assumption of polar symmetry 
requires us to use a reflective boundary condition on the 
symmetry axis, which causes material to slowly accumu- 
late in the polar regions. This makes it more difficult 
for the jet to blast through. However, in our earlier low 
resolution simulations, we have also simulated the pro- 
genitor evolution in the equatorial plane, thus avoiding 
the boundary condition problem. Although the details 
are underresolved, the jet still does not survive a massive 
WR shell. In this case the duration of the WR phase af- 
ter which the jet stalled was about 5000 yr, which may be 
an over-estimate due to the low resolution and therefore 
slower development of the shell. 

The duration of the WR phase cannot be determined 
exactly, because it is dependent on the development of 
the instabilities. This could explain the difference in the 
observed morphology of the two jet regions in Cas A 
(Fig-QJ- The northeast jet may have encountered a less 
massive portion of the shell, whereas the southwest jet 
may have been blocked by a density enhancement in the 
shell, thus creating a less distinct jet with a more dis- 
persed morphology. In fact, it is found that the velocity 
of the shock i n the west of Cas A is indeed lower than in 
the east (e.g. IVink et alj (|1998D ). 

The presence of a circumstellar shell also influences 



1000 2000 3000 4000 

Duration WR phase (yr) 

Fig. 5. — Energy required to accelerate the shell in the jet- 
region to velocities of 3iij e t/4, compared to the jet-energy (horizon- 
tal dashed line), for different implementations of the progenitor's 
mass loss history. Curve 'b' represents the simulation of the WR 
shell that was plotted in Figure [2] and used for the subsequent su- 
pernova simulations as plotted in Figure [4] Curve 'c' represents a 
simulation of the WR shell without initial density perturbations. 
Curve 'a' represents a simulation of WR shell with different initial 
density perturbations from the solid-line simulation. The dotted 
part of the line indicates where the jet starts to be effectively dis- 
persed into the rest of the remnant. Starting from a WR phase 
of about 1000-3000 years a jet with this energy does not protrude 
anymore. The vertical patch of line marks the time where analyt- 
ically the mass in the shell should be high enough to match the 
energy in the jet. 

the properties of the blast wave and reverse shock. For a 
longer WR phase the equatorial blast wave reaches fur- 
ther while the radius of the reverse shock decreases. This 
is summarized in Table Q] for an explosion with jet mass 
and energy of 3.8 x 10 29 g respectively 5.5 x 10 48 erg. Be- 
cause of the uncertainty in the observational constraints 
on the properties of the reverse shock, for all the scenar- 
ios in Table [T] the properties of the forward and reverse 
shock are within reasonable boundaries. It seems that 
the presence of a non-relativistic jet constrains the max- 
imal duration of the WR phase more strictly. 
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TABLE 1 



' WR SHELL AND FORWARD AND REVERSE SHOCK AND THE BLAST WAVE OF THE JET AFTER t = 








EjET = 5.5 X 10 48 ERG. 








WR Phase 


^shell 


LVJl W CLL K_l 


^forward 


^reverse 


"reverse 


^jet 


v \et 


(yr) 


(PC) 


(PC) 


(cm s _1 ) 


(pc) 


(cm s _1 ) 


(PC) 


(cm s -1 ) 








2.43 


5.43 X 10 s 


1.63 


3.16 x 10 8 


> 6 


1.40 X 10 M 


1569 


0.26 


2.44 


5.34 x 10 s 


1.28 


2.57 x 10 s 


4.25 


1.03 x 10 9 


1743 


0.28 


2.45 


5.33 x 10 s 


1.24 


2.27 x 10 8 


5.54 


1.24 x 10 9 


2091 


0.33 


2.49 


5.53 x 10 s 


1.07 


2.45 x 10 s 


> 6 


1.45 x 10 9 


2266 


0.35 


2.50 


5.93 x 10 s 


0.93 


9.90 x 10 7 


4.57 


9.90 x 10 s 


2440 


0.37 


2.47 


4.94 x 10 s 


0.90 


1.98 x 10 s 


3.58 


7.93 x 10 s 


2614 


0.40 


2.47 


5.24 x 10 8 


0.96 


1.88 x 10 8 


3.17 


7.56 x 10 s 


3486 


0.50 


2.49 


5.33 x 10 s 


0.95 


1.28 x 10 s 


2.73 


5.73 x 10 s 


From observations: 


2.5 


5.0 x 10 8 


1.4 


3.0 x 10 s 


> 3.8 





330 YR AND FOR 



Note. — Properties of the SNR and the jet depend on the CSM the blastwave runs into. The extent of the WR shell before the explosion 
is given in the second column. The amount of mass swept up in the shell is essentially the mass of the RSG wind within this radius. The 
WR wind does not contribute much mass. The bottom row shows the preferred value of the SNR parameter, from observational constraints. 
The jet puts tighter restrictions on the duration of the WR phase than the properties of the blast wave and reverse shock alone. 



4. CONCLUSIONS AND DISCUSSION 

We have simulated the evolution of axisymmetric 
ejecta, such as may result from a bipolar supernova ex- 
plosion, in the context of Cas A and the presence of a 
jet / counter-jet in this SNR. For the initial conditions 
we used a realistic progenitor evolution for a ~ 20 M Q 
star, consisting of a RSG wind, followed by a WR wind. 
We find that the presence of a WR shell limits the sur- 
vival of the jets. The survival depends critically on the 
energy available in the jet region and on the mass con- 
tained in the WR shell. The latter is determined by the 
duration of the WR phase and the properties of the pro- 
genitor winds. For the parameters chosen for the Cas A 
progenitor, we find that if the WR phase is longer than 
2000—5000 yr, the jets do not protrude through the shell, 
in which case the situation does not correspond to the 
presence of jets in Cas A. Therefore, either the progen- 
itor went through a very short WR phase, or it did not 
have one and exploded as a RSG. In general however, 
this means that, also if a SNR appears symmetric, the 
explosion may still have been accompanied by jets. 

In order to match the length of the observed jets of 
Cas A, an energy of at least 2.0 x 10 48 erg per jet is re- 
quired. In case a WR shell is formed, higher energies are 
needed. However, for a WR duration in excess of maxi- 
mally 5000 yr, the properties of the remnant and the jet 
do not match the observations. The upper limit we find 
corresponds to the findings of van Veelen & Langer (2008, 
in preparation), where they find that the properties of the 
forward and reverse shock do not agree with observations 
for a WR phase that lasts more than 5000 years. 

The question now arises if a scenario involving a very 
short WR phase is realistic. If the progenitor was not in 
a binary, the chance of having a very short WR phase is 
very low. The main reason for invoking a WR phase at all 
is that it explains the dumpiness in the remnant, the lack 
of hydrogen in the ejecta, the presence of metal clumps 
far out in the ejecta, and the high N/H ratio in the CSM. 
A single star at th e end of the RSG phase w ith a clumpy 
mass loss history ([Chevalier fc Oishill2003f) may also be 
able to partly explain the above, but requir e s eq ually 
coincidental circumstances. lLaming fc Hwangl (|2007h fa- 
vor a short WR phase because, based on the temperature 
and ionization age of the X-ray emitting gas, they find 
that the ejecta expanded in a bubble of around ~ 0.2 pc, 
in our case corresponding to a WR phase of about 1000 



yrs. They however do not take into account the presence 
of a shell around the bubble. 

As an alternative to the single star model, a model with 
a bin ary companion has been proposed (jYoung et al.l 
200(1). A common envelope (CE) phase in a close binary 
solves a number of problems: It explains the low ejecta 
mass in conjunction with a MS mass of ~ 20 Mq, and 
provides a natural explanation for a very short WR phase 
of the primary star (Podsiadlowski, private communica- 
tion). However, the details of CE evolution are not well 
understood and no companion star has currently been 
found. 

The simulated jets resemble the observed jet of Cas A 
in opening angle for a variety of parameters. This rein- 
forces the idea that the explosion itself was intrinsically 
bipolar. It remains an interesting question what mech- 
anism is responsible for such an asymmetric explosion, 
and whether it is related to other bipolar explosion phe- 
nomena such as X-ray flashes or LGRBs. Although our 
case does not resemble the relativistic sce narios as ar e 
invoked in models for gamma ray bursts (|Piranl 120051) . 
it does not seem unreasonable that once again rotation 
is involved in creating the asymmetry in a low-energy 
explosion like this. Unless the rotation was creat ed dur- 
ing the explosion (Blo ndin fc Mezzacappal l2~007t ). rota- 
tion may have also left an imprint on the CSM and a 
combination of asymmetric CSM and an asymmetric ex- 
plosion could have been responsible for the Cas A mor- 
phology. The observed point source near the center of 
the remnant is likely to be a neutron star. However, the 
absence of a bright pulsar wind nebula suggests that the 
present rotation peri od of the neutron star is relativel y 
low (> 160-330 ms lSeward fc Wandfl988L IVml[2T)0l . 
This seems at odds with a rapid rotation of the stellar 
core as a mechanism to create a bipolar explosion. How- 
ever, this discrepancy may be solved if the point source 
is a magnetar that has considerably slowed down. 

In conclusion, we would like to emphasize that the 
presence of jets in Cas A, together with considerable 
knowledge about nucleosynthesis yields, explosion en- 
ergy, and compact object, makes this SNR a unique ob- 
ject to investigate the mechanism behind bipolar explo- 
sions, and the type of progenitors that produce bipolar 
explosions. 
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APPENDIX 

LENGTH-ENERGY RELATION FOR A CONICAL AND CYLINDRICAL JET 

The length-energy relation for a jet has been derived in the approximation of a conical and cylindrical morphology 
for the jet. We use assumptions similar to those that lead to the well-known Sedov- Taylor solution for a spherical blast 
wave, i.e. a jet interior that is significantly over-pressured with respect to the surrounding medium and no significant 
radiative losses. For a jet with a constant opening angle, energy conservation gives: 



E = -Mv 2 + - . (Al) 

2 4tt 3 7 - 1 y ' 

We assume pressure equilibrium between the SNR interior and the shell, and take for the velocity the post-shock 
velocity, thus getting: 

(A2) 

1 /2 

Rl ' 2dR= {^lAj dt ' (A3) 

r <^a)" 3 (A4) 

/ Et 2\ 1/3 

RKlA5 {efA) ' (A5) 

Here 9 } is the jet's opening angle, E the jet-energy, R the length of the jet at time t, r the radius of the cone, p = A/R 2 
with A a normalization constant, and C 7 = (47)/((7 — 1)(7 + l) 2 ) = 45/32 for a polytropic index of 7 = 5/3. 

For a cylinder we use the same approach, but because of the different volume we get a different relation between the 
length and energy: 

E = C 7 m(^)\ (A6) 



dt J 



The mass of a cylinder-like jet is: 



M 

'Ro 

f Rl 2 A 



= f 1 nr 2 p dR (A7) 

J Rn 

nr 2 4^dR (A8) 



Energy conservation is then given by: 



,R-R fdR\ 2 E 
R — Rq ( dR^ 



B = C,-W— -2 — 7^3 = < (A10) 



sir {.*) ■ (A11) 

Substitution of x = with changing the boundaries from (R ,R) accordingly to (1, -^) leads to: 

x — 1 „ f dx\ 2 



.l — 1 dx , . „ „. 

^=\I—Tf (A13) 
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This can be reduced to quadrature by substituting x = y 2 : 



ry/R/Ra i 

J \/y 2 -^dy=-^/lt. 



The left side is a standard integral: 



r^/R/Ra ^ i . . 

J Vy 2 - ld v = -jvvv 2 - 1 - 2 ln [y + v y 2 - 1 J 



\J R/ Ra 



Thus the solution becomes: 

yVv 2 - 1 - in (y + Vy 2 - i) 

Restoring the original variables gives: 



— /— -1-lnl /— + /— -1 

-Ro V -Ro \ V -Ro V -Ro 



Ro 



For R^> Rq this can be approximated by: 



R = 




Et 2 R I I R . 



(A14) 



(A15) 



(A16) 



(A17) 



(A18) 



(A19) 



The second term on the right-hand side has only a very weak dependence on R and in the limit R ^> Ro it is 
negligable compared to the first term on the right-hand side. 
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